High-resolution KIR and HLA genotyping in three Chinese ethnic minorities reveals distinct origins.

Polymorphism of killer-cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands impacts the effector activity of cytotoxic NK cell and T cell subsets. Therefore, understanding the extent and implications of KIR and HLA class I genetic polymorphism across various populations is important for immunological and medical research. In this study, we conducted a high-resolution investigation of KIR and HLA class I diversity in three distinct Chinese ethnic minority populations. We studied the She, Yugur, and Tajik, and compared them with the Zhejiang Han population (Zhe), which represents the majority Southern Han ethnicity. Our findings revealed that the Tajik population exhibited the most diverse KIR copy number, allele, and haplotype diversity among the four populations. This diversity aligns with their proposed ancestral origin, closely resembling that of Iranian populations, with a relatively higher presence of KIR-B genes, alleles, and haplotypes compared with the other Chinese populations. The Yugur population displayed KIR distributions similar to those of the Tibetans and Southeast Asians, whereas the She population resembled the Zhe and other East Asians, as confirmed by genetic distance analysis of KIR. Additionally, we identified 12.9% of individuals across the three minority populations as having KIR haplotypes characterized by specific gene block insertions or deletions. Genetic analysis based on HLA alleles yielded consistent results, even though there were extensive variations in HLA alleles. The observed variations in KIR interactions, such as higher numbers of 2DL1-C2 interactions in Tajik and Yugur populations and of 2DL3-C1 interactions in the She population, are likely shaped by demographic and evolutionary mechanisms specific to their local environments. Overall, our findings offer valuable insights into the distribution of KIR and HLA diversity among three distinct Chinese ethnic minority populations, which can inform future clinical and population studies.

Donor KIR2DL1 Allelic Polymorphism Influences Posthematopoietic Progenitor Cell Transplantation Outcomes in the T Cell Depleted and Reduced Intensity Conditioning Setting.

The majority of established KIR clinical assessment algorithms used for donor selection for hematopoietic progenitor cell transplantation (HPCT) evaluate gene content (presence/absence) of the KIR gene complex. In comparison, relatively little is known about the impact of KIR allelic polymorphism. By analyzing donors of T cell depleted (TcD) reduced intensity conditioning (RIC) HPCT, this study investigated the influence on post-transplant outcome of 2 polymorphic residues of the inhibitory KIR2DL1. The aim of this study was to expand upon existing research into the influence of KIR2DL1 allelic polymorphism upon post-transplant outcome. The effects of allele groups upon transplant outcomes were investigated within a patient cohort using a defined treatment protocol of RIC with TcD. Using phylogenetic data, KIR2DL1 allelic polymorphism was categorized into groups on the basis of variation within codons 114 and 245 (positive or negative for the following groups: KIR2DL1*002/001g, KIR2DL1*003, KIR2DL1*004g) and the identification of null alleles. The influence of these KIR2DL1 allele groups in hematopoietic progenitor cell transplantation (HPCT) donors was assessed in the post-transplant data of 86 acute myelogenous leukemia patients receiving RIC TcD HPCT at a single center. KIR2DL1 allele groups in the donor significantly impacted upon 5-year post-transplant outcomes in RIC TcD HPCT. Donor KIR2DL1*003 presented the greatest influence upon post-transplant outcomes, with KIR2DL1*003 positive donors severely reducing 5-year post-transplant overall survival (OS) compared to those receiving a transplant from a KIR2DL1*003 negative donor (KIR2DL1*003 pos versus neg: 27.0% versus 60.0%, P = .008, pc = 0.024) and disease-free survival (DFS) (KIR2DL1*003 pos versus neg: 23.5% versus 60.0%, P = .004, pc = 0.012), and increasing 5-year relapse incidence (KIR2DL1*003 pos versus neg: 63.9% versus 27.2%, P = .009, pc = 0.027). KIR2DL1*003 homozygous and KIR2DL1*003 heterozygous grafts did not present significantly different post-transplant outcomes. Donors possessing the KIR2DL1*002/001 allele group were found to significantly improve post-transplant outcomes, with donors positive for the KIR2DL1*004 allele group presenting a trend towards improvement. KIR2DL1*002/001 allele group (KIR2DL1*002/001g) positive donors improved 5-year OS (KIR2DL1*002/001g pos versus neg: 56.4% versus 27.2%, P = .009, pc = 0.024) and DFS (KIR2DL1*002/001g pos versus neg: 53.8% versus 25.5%, P = .018, pc = 0.036). KIR2DL1*004 allele group (KIR2DL1*004g) positive donors trended towards improving 5-year OS (KIR2DL1*004g pos versus neg: 53.3% versus 35.5%, P = .097, pc = 0.097) and DFS (KIR2DL1*004g pos versus neg: 50.0% versus 33.9%, P = .121, pc = 0.121), and reducing relapse incidence (KIR2DL1*004g pos versus neg: 33.1% versus 54.0%, P = .079, pc = 0.152). The presented findings suggest donor selection algorithms for TcD RIC HPCT should consider avoiding KIR2DL1*003 positive donors, where possible, and contributes to the mounting evidence that KIR assessment in donor selection algorithms should reflect the conditioning regime protocol used.

Exceptional diversity of KIR and HLA class I in Egypt.

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases.

HLA-DPB1*13:01 associates with enhanced, and KIR2DS4*001 with diminished protection from developing severe COVID-19.

Extreme polymorphism of HLA and killer-cell immunoglobulin-like receptors (KIR) differentiates immune responses across individuals. Additional to T cell receptor interactions, subsets of HLA class I act as ligands for inhibitory and activating KIR, allowing natural killer (NK) cells to detect and kill infected cells. We investigated the impact of HLA and KIR polymorphism on the severity of COVID-19. High resolution HLA class I and II and KIR genotypes were determined from 403 non-hospitalized and 1575 hospitalized SARS-CoV-2 infected patients from Italy collected in 2020. We observed that possession of the activating KIR2DS4*001 allotype is associated with severe disease, requiring hospitalization (OR = 1.48, 95% CI 1.20-1.85, pc = 0.017), and this effect is greater in individuals homozygous for KIR2DS4*001 (OR = 3.74, 95% CI 1.75-9.29, pc = 0.003). We also observed the HLA class II allotype, HLA-DPB1*13:01 protects SARS-CoV-2 infected patients from severe disease (OR = 0.49, 95% CI 0.33-0.74, pc = 0.019). These association analyses were replicated using logistic regression with sex and age as covariates. Autoantibodies against IFN-α associated with COVID-19 severity were detected in 26% of 156 hospitalized patients tested. HLA-C*08:02 was more frequent in patients with IFN-α autoantibodies than those without, and KIR3DL1*01502 was only present in patients lacking IFN-α antibodies. These findings suggest that KIR and HLA polymorphism is integral in determining the clinical outcome following SARS-CoV-2 infection, by influencing the course both of innate and adaptive immunity.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

Unraveling the Phenotypic States of Human innate-like T Cells: Comparative Insights with Conventional T Cells and Mouse Models.

The "innate-like" T cell compartment, known as Tinn, represents a diverse group of T cells that straddle the boundary between innate and adaptive immunity, having the ability to mount rapid responses following activation. In mice, this ability is acquired during thymic development. We explored the transcriptional landscape of Tinn compared to conventional T cells (Tconv) in the human thymus and blood using single cell RNA sequencing and flow cytometry. We reveal that in human blood, the majority of Tinn cells, including iNKT, MAIT, and Vδ2+Vγ9+ T cells, share an effector program characterized by the expression of unique chemokine and cytokine receptors, and cytotoxic molecules. This program is driven by specific transcription factors, distinct from those governing Tconv cells. Conversely, only a fraction of thymic Tinn cells displays an effector phenotype, while others share transcriptional features with developing Tconv cells, indicating potential divergent developmental pathways. Unlike the mouse, human Tinn cells do not differentiate into multiple effector subsets but develop a mixed type I/type III effector potential. To conduct a comprehensive cross-species analysis, we constructed a murine Tinn developmental atlas and uncovered additional species-specific distinctions, including the absence of type II Tinn cells in humans, which implies distinct immune regulatory mechanisms across species. The study provides insights into the development and functionality of Tinn cells, emphasizing their role in immune responses and their potential as targets for therapeutic interventions.

Allelic variation of KIR and HLA tunes the cytolytic payload and determines functional hierarchy of NK cell repertoires.

The functionality of natural killer (NK) cells is tuned during education and is associated with remodeling of the lysosomal compartment. We hypothesized that genetic variation in killer cell immunoglobulin-like receptor (KIR) and HLA, which is known to influence the functional strength of NK cells, fine-tunes the payload of effector molecules stored in secretory lysosomes. To address this possibility, we performed a high-resolution analysis of KIR and HLA class I genes in 365 blood donors and linked genotypes to granzyme B loading and functional phenotypes. We found that granzyme B levels varied across individuals but were stable over time in each individual and genetically determined by allelic variation in HLA class I genes. A broad mapping of surface receptors and lysosomal effector molecules revealed that DNAM-1 and granzyme B levels served as robust metric of the functional state in NK cells. Variation in granzyme B levels at rest was tightly linked to the lytic hit and downstream killing of major histocompatibility complex-deficient target cells. Together, these data provide insights into how variation in genetically hardwired receptor pairs tunes the releasable granzyme B pool in NK cells, resulting in predictable hierarchies in global NK cell function.

Polymorphic KIR3DL3 expression modulates tissue-resident and innate-like T cells.

Most human killer cell immunoglobulin-like receptors (KIR) are expressed by natural killer (NK) cells and recognize HLA class I molecules as ligands. KIR3DL3 is a conserved but polymorphic inhibitory KIR recognizing a B7 family ligand, HHLA2, and is implicated for immune checkpoint targeting. The expression profile and biological function of KIR3DL3 have been somewhat elusive, so we searched extensively for KIR3DL3 transcripts, revealing highly enriched expression in γδ and CD8+ T cells rather than NK cells. These KIR3DL3-expressing cells are rare in the blood and thymus but more common in the lungs and digestive tract. High-resolution flow cytometry and single-cell transcriptomics showed that peripheral blood KIR3DL3+ T cells have an activated transitional memory phenotype and are hypofunctional. The T cell receptor (TCR) usage is biased toward genes from early rearranged TCR-α variable segments or Vδ1 chains. In addition, we show that TCR-mediated stimulation can be inhibited through KIR3DL3 ligation. Whereas we detected no impact of KIR3DL3 polymorphism on ligand binding, variants in the proximal promoter and at residue 86 can reduce expression. Together, we demonstrate that KIR3DL3 is up-regulated alongside unconventional T cell stimulation and that individuals may vary in their ability to express KIR3DL3. These results have implications for the personalized targeting of KIR3DL3/HHLA2 checkpoint inhibition.

KIR-HLA interactions extend human CD8+ T cell lifespan in vivo.

BACKGROUNDThere is increasing evidence, in transgenic mice and in vitro, that inhibitory killer cell immunoglobulin-like receptors (iKIRs) can modulate T cell responses. Furthermore, we have previously shown that iKIRs are an important determinant of T cell-mediated control of chronic viral infection and that these results are consistent with an increase in the CD8+ T cell lifespan due to iKIR-ligand interactions. Here, we tested this prediction and investigated whether iKIRs affect T cell lifespan in humans in vivo.METHODSWe used stable isotope labeling with deuterated water to quantify memory CD8+ T cell survival in healthy individuals and patients with chronic viral infections.RESULTSWe showed that an individual's iKIR-ligand genotype was a significant determinant of CD8+ T cell lifespan: in individuals with 2 iKIR-ligand gene pairs, memory CD8+ T cells survived, on average, for 125 days; in individuals with 4 iKIR-ligand gene pairs, the memory CD8+ T cell lifespan doubled to 250 days. Additionally, we showed that this survival advantage was independent of iKIR expression by the T cell of interest and, further, that the iKIR-ligand genotype altered the CD8+ and CD4+ T cell immune aging phenotype.CONCLUSIONSTogether, these data reveal an unexpectedly large effect of iKIR genotype on T cell survival.FUNDINGWellcome Trust; Medical Research Council; EU Horizon 2020; EU FP7; Leukemia and Lymphoma Research; National Institute of Health Research (NIHR) Imperial Biomedical Research Centre; Imperial College Research Fellowship; National Institutes of Health; Jefferiss Trust.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures.

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

The impact of HLA polymorphism on herpesvirus infection and disease.

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

Host KIR/HLA-C Genotypes Determine HIV-Mediated Changes of the NK Cell Repertoire and Are Associated With Vpu Sequence Variations Impacting Downmodulation of HLA-C.

NK cells play a pivotal role in viral immunity, utilizing a large array of activating and inhibitory receptors to identify and eliminate virus-infected cells. Killer-cell immunoglobulin-like receptors (KIRs) represent a highly polymorphic receptor family, regulating NK cell activity and determining the ability to recognize target cells. Human leukocyte antigen (HLA) class I molecules serve as the primary ligand for KIRs. Herein, HLA-C stands out as being the dominant ligand for the majority of KIRs. Accumulating evidence indicated that interactions between HLA-C and its inhibitory KIR2DL receptors (KIR2DL1/L2/L3) can drive HIV-1-mediated immune evasion and thus may contribute to the intrinsic control of HIV-1 infection. Of particular interest in this context is the recent observation that HIV-1 is able to adapt to host HLA-C genotypes through Vpu-mediated downmodulation of HLA-C. However, our understanding of the complex interplay between KIR/HLA immunogenetics, NK cell-mediated immune pressure and HIV-1 immune escape is still limited. Therefore, we investigated the impact of specific KIR/HLA-C combinations on the NK cell receptor repertoire and HIV-1 Vpu protein sequence variations of 122 viremic, untreated HIV-1+ individuals. Compared to 60 HIV-1- controls, HIV-1 infection was associated with significant changes within the NK cell receptor repertoire, including reduced percentages of NK cells expressing NKG2A, CD8, and KIR2DS4. In contrast, the NKG2C+ and KIR3DL2+ NK cell sub-populations from HIV-1+ individuals was enlarged compared to HIV-1- controls. Stratification along KIR/HLA-C genotypes revealed a genotype-dependent expansion of KIR2DL1+ NK cells that was ultimately associated with increased binding affinities between KIR2DL1 and HLA-C allotypes. Lastly, our data hinted to a preferential selection of Vpu sequence variants that were associated with HLA-C downmodulation in individuals with high KIR2DL/HLA-C binding affinities. Altogether, our study provides evidence that HIV-1-associated changes in the KIR repertoire of NK cells are to some extent predetermined by host KIR2DL/HLA-C genotypes. Furthermore, analysis of Vpu sequence polymorphisms indicates that differential KIR2DL/HLA-C binding affinities may serve as an additional mechanism how host genetics impact immune evasion by HIV-1.

Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health.

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Human herpesvirus diversity is altered in HLA class I binding peptides.

Herpesviruses are ubiquitous, genetically diverse DNA viruses, with long-term presence in humans associated with infrequent but significant pathology. Human leukocyte antigen (HLA) class I presents intracellularly derived peptide fragments from infected tissue cells to CD8+ T and natural killer cells, thereby directing antiviral immunity. Allotypes of highly polymorphic HLA class I are distinguished by their peptide binding repertoires. Because this HLA class I variation is a major determinant of herpesvirus disease, we examined if sequence diversity of virus proteins reflects evasion of HLA presentation. Using population genomic data from Epstein­Barr virus (EBV), human cytomegalovirus (HCMV), and Varicella­Zoster virus, we tested whether diversity differed between the regions of herpesvirus proteins that can be recognized, or not, by HLA class I. Herpesviruses exhibit lytic and latent infection stages, with the latter better enabling immune evasion. Whereas HLA binding peptides of lytic proteins are conserved, we found that EBV and HCMV proteins expressed during latency have increased peptide sequence diversity. Similarly, latent, but not lytic, herpesvirus proteins have greater population structure in HLA binding than nonbinding peptides. Finally, we found patterns consistent with EBV adaption to the local HLA environment, with less efficient recognition of EBV isolates by high-frequency HLA class I allotypes. Here, the frequency of CD8+ T cell epitopes inversely correlated with the frequency of HLA class I recognition. Previous analyses have shown that pathogen-mediated natural selection maintains exceptional polymorphism in HLA residues that determine peptide recognition. Here, we show that HLA class I peptide recognition impacts diversity of globally widespread pathogens.

Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1.

Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.

The polymorphism L412F in TLR3 inhibits autophagy and is a marker of severe COVID-19 in males.

The polymorphism L412F in TLR3 has been associated with several infectious diseases. However, the mechanism underlying this association is still unexplored. Here, we show that the L412F polymorphism in TLR3 is a marker of severity in COVID-19. This association increases in the sub-cohort of males. Impaired macroautophagy/autophagy and reduced TNF/TNFα production was demonstrated in HEK293 cells transfected with TLR3L412F-encoding plasmid and stimulated with specific agonist poly(I:C). A statistically significant reduced survival at 28 days was shown in L412F COVID-19 patients treated with the autophagy-inhibitor hydroxychloroquine (p = 0.038). An increased frequency of autoimmune disorders such as co-morbidity was found in L412F COVID-19 males with specific class II HLA haplotypes prone to autoantigen presentation. Our analyses indicate that L412F polymorphism makes males at risk of severe COVID-19 and provides a rationale for reinterpreting clinical trials considering autophagy pathways.Abbreviations: AP: autophagosome; AUC: area under the curve; BafA1: bafilomycin A1; COVID-19: coronavirus disease-2019; HCQ: hydroxychloroquine; RAP: rapamycin; ROC: receiver operating characteristic; SARS-CoV-2: severe acute respiratory syndrome coronavirus 2; TLR: toll like receptor; TNF/TNF-α: tumor necrosis factor.

Estimating HLA haplotype frequencies from homozygous individuals - A Technical Report.

We estimated HLA haplotype frequencies based on individuals homozygous for 4, 5 or 6 loci. Validation of our approach using a sample of over 3.4 million German individuals was successful. Compared to an expectation-maximization algorithm, the errors were larger. However, our approach allows the unequivocal detection of rare haplotypes.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING.

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.